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FATIGUE  LITE  PREDICTION 
FOR  SIMULTANEOUS 
STRESS  AND  STRENGTH  VARIANCES 

Simple  closed  form  expressions  have  been  found  Co  accurately 
predict  the  fatigue  life  of  structures  subjected  to  sinusoidal  or  random 
stresses  where  Che  applied  stress  and  the  material's  strength  are 
simultaneous  random  variables.  With  appropriate  parameter  value  changes 
the  same  equations  accurately  apply  to  both  the  low  cycle  (inelastic)  and 
high  cycle  (elastic)  fatigue  regions.  These  equations  are  in  familiar 
engineering  terms.  Comparisons  between  analytical  predictions  and  empirical 
results  have  shown  Co  be  good  whenever  such  comparisons  were  made. 

INTRODUCTION 

Many  closed  form  analytical  expressions  have  previously  been  derived 
to  predict  structural  fatigue  life  and  mechanical  reliability  for  both 
sinusoidally  and  randomly  applied  stresses  and  strains  {1}  {2}  {3}  {4}  {5}. 
These  expressions  have  been  shown  to  be  simple,  practical  and  accurate.  They 
apply  to  single  and  multi-degeee-of-freedom  systems,  to  single  level  or  step- 
stress  load  situations,  and  to  both  low  and  high  cycle  fatigue  regions. 
Fracture  Mechanics  effects  are  included.  In  all  of  these  cases  the  stress/ 
strength  and  strain/ductility  parameters  were  treated  as  random  variables 
independently,  not  simultaneously. 

In  most  practical  cases  the  stress/strength  parameters  are  simultaneous 
random  variables.  Stresses  vary  from  part  to  part  and  subassembly  to  sub- 
assembly  due  to  dimensional  and  geometrical  differences  between  parts, 
fabrication  and  assembly  variances,  and  structural  damping  and  stiffness 
variances  of  adjacent  structures.  Strengths  vary  because  materials’  fatigue 
curves  are  a  scatterband  of  failure  points,  not  single  lines. 


1 


APPROACH  SUMMARY 


An  attempt  to  rigorously  derive  a  fatigue  life  expression  with  the 
stress/strength  parameters  treated  as  simultaneous  random  variables  was 
unsuccessful  in  that  the  final  expression  was  exceedingly  complex. 

Therefore  a  different  approach  was  evaluated.  This  approach  modified  the 
variable  strength  fatigue  life  expression  Cl}  by  adding  the  stress  (6)  and 
strength  (A)  standard  deviations  in  the  mean-square  sense  and  substituting 
the  resulting  standard  deviation  (  ¥  - \Jt T  +  6^  )  in  place  of  the  strength 
standard  deviation  term  (A) .  The  reasoning  behing  this  approach  was  as 
follows:  Fatigue  failure  occurs  when  stress  exceeds  strength  regardless  of 
whether  the  stress  is  "too  high"  or  the  strength  is  "too  low".  Both 
deviations  from  nominal  cause  a  reduction  in  fatigue  life.  Since  the  standard 
deviations  of  stress  and  strength  are  independent  of  each  other,  they  should 
be  added  in  the  mean-square  sense.  This  approach,  as  judged  by  Monte-Carlo 
simulation  techniques,  gives  somewhat  accurate  results  but  not  as  accurate 
as  hoped  for. 

Accuracy  was  improved  by  multiplying  the  stress  standard  deviation  (6) 

2/S 

by  the  term  (21^)  .  Njq  is  the  median  stress  cycles  to  failure.  It  is  the 

fatigue  life  if  the  analysis  is  done  deterministically  (i.e.  if  A  and  5  are 
zero) .  3  is  the  slope  parameter  of  the  materials  "S-N"  fatigue  curve.  This 

term  made  the  entire  expression  almost  identical  to  the  rigorously  derived 
equation  for  the  case  of  A*  0. 

Accuracy  was  further  improved  in  the  region  of  early  fatigue  failures 

1/3 

by  subtracting  the  term  (21^)  A<3. 

<\J  28  -  tt/S  ' 
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APPROACH  SUMMARY  (Cont'd) 


The  portion^/28  —  ir/6 '  was  required  to  provide  accuracy  for  both  the  low 
and  high  cycle  fatigue  regions  and  for  brittle  and  ductile  materials.  This 
worsened  the  accuracy  in  the  region  of  the  late  failures.  The  above  term 
needed  to  be  added  instead  of  subtracted  in  that  region  (i.e.  a  sign  change 
for  N  >  Njj) .  This  worsened  the  accuracy  in  the  middle  failure  region.  The 
multiplying  term  5*2  erf 
regions. 


"20(i  •  j 


restored  accuracy  to  all  failure 


The  resultant  standard  deviation  term  is 


/A2  +  (21^) 2/6  62  +  §  (2Nm)1,/0A,s 


Accuracy  of  the  above  expressions  was  judged  by  comparison  to  Monte- 
Carlo  simulation  results.  The  Monte  Carlo  simulation  technique  had  its 
accuracy  and  practicality  checked  by  comparing  its  results  with  those  known 
to  be  theoretically  correct  and  with  available  empirical  results. 


Fatigue  life  is  expressed  in  terms  of  probability  of  failure  as  a 
function  of  applied  stress  cycles  and  both  average  and  minimum  cycles  to 
first  failure.  For  the  most  part  data  is  presented  in  the  form  of  histo¬ 
grams  of  cycles  to  failure  because  of  the  histogram's  sensitivity  to 


differences  between  theoretical  and  tallied  results. 


SUMMARY  OF  RESULTS 


The  single  expression  for  ¥  provides  accurate  fatigue  life  results  for 
ductile  and  brittle  materials,  over  the  early  and  late  failure  regions 
ofr  both  low  (inelastic)  and  high  (elastic)  cycle  fatigue  situations.  All 
of  the  fatigue  life  and  mechanical  reliability  equations  in  references 
through  {5}  that  originally  applied  to  cases  where  strength  alone  was  the 
random  variable,  can  be  used  for  simultaneous  stress/strength  variances  by 
substituting  ¥  for  A. 

The  Monte  Carlo  simulation  technique  was  judged  to  be  both  accurate  and 
practical  due  to  good  comparisons  with  results  known  to  be  theoretically 
correct  and  with  available  empirical  results. 


FATIGUE  CURVE  REPRESENTATION 


Modern  fatigue  curve  representation  is  as  shown  in  figures  1  and  2  {6}  {7} 
Figure  1  is  a  plot  of  "true"  stress  amplitude  versus  reversals  to  failure 
for  1020  HR  steel.  It  covers  both  the  low  cycle  (plastic  or  inelastic) 
and  the  high  cycle  (elastic)  fatigue  regions.  The  fatigue  curve  is  a 
single  straight  line.  "True"  stress  is  defined  as  the  applied  load 
divided  by  the  actual  cross-section  area,  which  becomes  less  than  the 
original  area  as  the  load  is  increased.  "True"  stress  is  contrasted  with 
"engineering"  stress  which  is  defined  as  the  applied  load  divided  by  the 
original  cross-section  area.  '  Life  is  in  terms  of  reversals  to  failure  or 
twice  the  cycles  to  failures  Nf;  there  being  two  reversals  for  each  stress 
cycle.  The  fatigue  strength  coefficient  cr^  can  be  thought  of  as  being 
approximately  equal  to  the  "true"  ultimate  strength  of  the  material.  The 
fatigue  strength  exponent  can  be  thought  of  as  a  slope  parameter. 

Figure  2  is  a  plot  of  the  same  failure  data  as  figure  1  except  the  ordinate 
is  expressed  as  "true"  strain  amplitude.  The  strain-life  curve  is  the  sum 
of  the  plastic  and  elastic  strain-life  curves.  E  is  the  modulus  of 
elasticity.  is  the  fatigue  ductility  coefficient.  It  can  be  thought  of 

as  a  measure  of  the  material's  ductility.  The  fatigue  ductility  exponent  c 
has  a  value  of  approximately  -0.5  for  most  structural  materials.  The  fatigue 
strength  exponent  b  takes  on  values  of  approximately  -0.1  for  ductile 
materials  to  -0.05  for  brittle  materials.  Fatigue  curve  data  for  many 
materials  is  found  in  reference  {7}  .  The  strain-life  curve  of  figure  2 
shows  that  the  plastic  strain-life  predominates  below  approximately  10^ 
cycles;  whereas  the  elastic  strain-life  curve  predominates  above  10^  cycles. 
The  transition  cycles  varies  widely  for  different  materials. 
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Figure  3  is  the  plastic  strain-life  curve  of  figure  2.  Figure  4  is  the 
elastic  strain- life  curve  of  figure  2.  For  the  analyses  of  this  paper  the 
strain-life  curves  similar  to  those  of  figures  3  and  4  will  be  used  for 
various  materials.  The  general  form  to  be  used  is  shown  in  Figure  5.  In 
the  plastic  region  6  *  2  for  most  structural  materials,  eu  *  tf ^  and  e  *  Aep/2 
(compare  with  figure  3) .  In  the  elastic  region  6  a  8  for  ductile  materials  and 
8s  20  for  very  brittle  materials,  ey  ■  a  f/E  and  e  ■  ^Ee/2  (compare  with 
figure  4.)  The  curve  of  Figure  5  is  expressed  as 


Nf 


(1) 


where  *  cycles  to  failure 

e  *  applied  strain  amplitude,  inches/ inch 
ep  *  "y- intercept",  in/ in 

£,,  represents  the  material's  ductility  in  the  plastic  region  and  the  material's 
strength  in  the  elastic  region.  Equation  (1)  and  the  single  line  of  figure  5 
represent  a  deterministic  fatigue  curve.  Actual  fatigue  curves  are  scatter- 
bands  of  failure  points.  The  single  line  represents  the  median.  The  scatter- 
band  of  points  can  be  represented  by  letting  eu  in  equation  (1)  become  a 
Gaussian  random  variable  with  mean  value  eu  and  standard  deviation  Ae.  A 
random  variable  applied  strain  amplitude  can  also  be  represented  as  a  Gaussian 
random  variable  with  mean  value  e  and  standard  deviation  6  .  Equation  (1) 
then  becomes 

(2) 


where  N_ 


median  cycles  to  failure 
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average  value  of  e 
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Equation  (2)  can  equivalently  be  expressed  as  follows: 
In  the  low  cycle  region: 


In  the  high  cycle  region: 


0 


.  _  '/A 

>  I;=A2 


Ao 

"IF 


where  Acr_  applied  stress  amplitude  .  as  I 
2  7 


(a't  \ 

0 

N  - 

m 

For  a  random  applied  stress  of  tms 

r  rn 

"-1/0-/ 

2  °  f 

1 

V*  * 

\/~r 

<*)  J 

1/0 


N 


m 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 
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The  standard  deviations  associated  with  strain-life  fatigue  curves  are  A 


and  tS£  for  the  material's  ductility  and  applied  strain  amplitude  respectively. 


These  are  in  3train  units  of  inches  per  inch.  The  corresponding  standard 
deviations  associated  with  stress-life  fatigue  curves  are  as  follows: 


A  »  E  A 

z 

(9) 

0 

5  »  E  5 

(10) 

/v 

where  E  »  modulus  of  elasticity,  psi 
A,  5  v  psi 
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ANALYTICAL  DERIVATION 

The  derivation  of  the  fatigue  life  expressions  begins  with  the  derivations 
of  equations  for  the  probability  density  function  of  cycles  to  failures  p(Nj) 
and  the  probability  of  failure  at  N  applied  stress  cycles  F(N).  Appendix  A 
describes  the  derivation  of  p(Nf)  simultaneous  variations  in  eu  and  e.  It  is 


p(Nf) 


1/S  1/S  “1 

2  Nf 


8  Ae  6e  ^ 


t  VT 


-(h  -tv) 
r 


2  e  erf  (cxj) 


-h2/r 


+  TT  e 


(ID 


2h  erf  (f^)  ' 

J\f? 

where  the  variables  a^,  h,  r  and  v  are  complicated  functions  of  N^. 


F(N)  *  Probability  that  Ng  >  N 


F(N) 


P  (Nf)dNf 


(12) 


It  can  be  seen  by  examining  equations  (11)  and  (12)  that  finding  a  simple 
closed  form  expression  for  F(N)  does  not  appear  likely.  Without  a  simple 
expression  for  F(N)  the  derivation  of  the  average  cycles  to  first  failure 


N.  and  the  minimum  cycles  to  first  failure  cannot  proceed. 

MIN 
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SIMULATION  TECHNIQUE 


A  Monte  Carlo  technique  was  used  as  the  simulation  method  for  judging 

the  accuracy  of  the  proposed  fatigue  life  expressions.  From  equation  (1) 

8 


N 

f 


(1) 


A  sample  of  the  random  variable  N^  is  generated  by  generating  a  sample  of 
ey  and  £$  then  performing  the  operation  indicated  by  equation  (1).  Each 
sample  of  eu  is  drawn  from  a  Gaussian  distribution  of  mean  value  Ey  and 
standard  deviation  Ae.  Each  sample  of  e  is  similarly  drawn  from  a  Gaussian 
distribution  of  mean  e  and  standard  deviation  5£ .  Negative  values  of  eu 
and  e  are  discarded.  The  samples  of  Nf  are  sorted  and  stored  in  array  bins 
according  to  the  sample's  value.  The  quantity  of  Nj  samples  that  fall  into 
each  bin  is  summed  and  scored.  A  printout  of  the  quantity  of  samples  in  each 
bin  of  the  array  represents  a  histogram  of  Nj  for  specific  values  of  Ag,  e, 
6£  and  B- 


COMPARISON  OF  SIMULATION  AND  THEORETICAL  RESULTS 


Theoretical  results  are  obtained  as  follows: 

Define  N8  -  e/  (13) 

N8  is  the  numerator  of  equation  (1).  Refer  to  equation  (A- 2)  in 
Appendix  A. 

N8  -  x 
Therefore 


P(N8) 


1/8  -1 

(N8) 


exp 


Figure  6  shows  a  graphic  illustration  of  the  eu  -  N8 
transformation.  This  illustrates  the  reason  for  the 


(14) 

mathematical 
N8  histogram  shape. 


F(N) 


F(N) 


N 


s 


p(N8)dN8 


1/8 

N 


0.5  +  erf 


(15) 

(16) 


A  histogram  array  bin  quantity  q  for  a  bin  that  extends  from  Na  to  is 
q  «{F(Nb)  -  F(Na)^  S  (17) 

where  S  ■  total  N8  sample  size 


Table  I  shows  the  values  of  the  parameters  for  several  N8  histograms.  A 
wide  range  of  8  values  was  chosen. 

Refer  to  equation  (14) .  Note  that  a  value  of  S  ■  1  should  give  a  Gaussian 

_  / 

histogram  (i.e.  Case  [3J  ).  Figures  7,  8  and  9  show  good  agreement  between 

equation  (17)  and  the  program  generated  histogram. 
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TABLE  I  N8  HISTOGRAM  DATA 


CASE 

EH 

Ae 

c 

3 

N8 

0.2 

0.3 

0.00079 

0.00632 

2 

0.041 

©' 

0.0185 

0.000925 

0.0001726 

0.0034515 

9.6 

2.5E-17 

©' 

0.0185 

0.000925 

0.004387 

0.00022 

12.1 

1.3E-21 

s' 

0.2 

0.03 

0 

0.0024 

1 

0.2 

>  A£  ,  5£  ,  e  ,  N8  ,  ^  IN/IN 
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Figures  10,  11  and  12  also  show  good  agreement  between  theoretical  and 
Monte  Carlo  results.  These  curves  also  show  the  expected  skewing  effect 
of  large  8  values. 

g 

Similar  results  would  be  obtained  for  generating  samples  for  e  ,  the 
demoniator  of  equation  (1) . 

In  later  sections  comparisons  will  be  made  between  Monte  Carlo  results  and 
theoretical  ones  for  where  £y  and  e  are  random  variables  independently. 

In  all  cases  it  will  be  seen  that  the  Monte  Carlo  simulation  technique  is 
accurate  compared  to  theoretical  equations. 
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PROPOSED  FATIGUE  LIFE  EXPRESSIONS 


Equations  (2)  through  (10)are  proposed  to  calculate  the  median  cycles  to 
failure  Nm  in  terms  of  either  stress /strength  or  strain/ductility  parameters. 
Equations  (9)  and  (10)  relate  the  standard  deviations  of  stress  and  strain. 
These  equations  are  to  be  used  in  the  following  proposed  fatigue  life 
expressions: 


F(N)  ■  Probability  of  failing  at  N  applied  stress  cycles 


erf  (a)  - 


-y  /  2 

e  dy 


/  ,  2/3  2  1/6 

-  /V  +  (2V  Se  +  §  (2V  Ae  5e  in/in 


26  -  tt/8 


§  -  2  erf 


[2°  ft.  -  ?] 


Figure  13  is  a  plot  of  §  versus  N/Nm. 
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The  above  equations  have  been  included  along  with  the  previously  described 
Monte  Carlo  simulation  technique  into  one  program  PL-1  which  is  written  in 
Basic  Language. 


PL-1 


SIMULATION  COMPUTER 
PROGRAM  LISTING 


10  REM  N1 =AVG  DUCTILITY ( IN/ IN) 

20  RE1T~D1*5TD'  DEV  DUCTTL I T YTIN/IN) 

30  REM  N2=AVG  APPLIED  STRAIN ( IN/ IN) 

40  REM"DZ=STD  DEV  AFPUI  ED~5T*AfN(  IN7Tn1 - 

50  REM 

Tii  REM- B"i  =BET A  ( SLOPE  PA^SmTTBRI 
70  REM  W=BIN  WIDTH  (CYCLES" 

W  REM~S=SAWPUF'S  IZE - 

90  REM  N8=GENERATED  CYCLES  AT  FAILURE 

“00-  REM - (-RANDOM- VARIABLE) - 

110  REM 

120  HtMX4= YAllied  aVG  dJCTITITY 
130  REM  D5=TALLIED  STD  DEV  DUCTILITY 

“40-  REM"T5=T  ALL  I  ED  "APPLIED  STRATA - 

150  REM  D6=TALLIED  STD  DEV  APPLIED  STRAIN 

T60  REM - 

170  REM  F1=FAILURE  PROBABILITY  USING  TALLIED  PARAMETERS 

130'  NTaTT35 - 

190  01= . 000925 


"200  02  g.  00022 - 

210  N2=. 004337 

220  PR INI - 

230  PRINT  "DESIRED  PARAMETERS:" 

Z4fi  PRINT  "DUCTILlTYrAVGf STD  DEV:" 

250  PRINT  Nl»Dl 

TS0  PRINT  "APPLIED  3  TR  A I NTAVG  TSTirDW 
270  PRINT  N2.D2 

-30"  PRTWT - 

290  L= 1E7 

00"N3="0 - 

310  33=0 

—  0  T3*0 - 

330  S9=0 


) 


40  T7=0 


350  S1=0 


'360  Kgl - 

370  Bl=9.6 

330  SgTTTTgn - 

390  B3=2*B2 

400  W»5E5 - 7 - 

410  3=10000 
420  A1=0. 2543296 
430  A2=-0. 2844967 

-440  -fr?*T.-4 2 1414 - 

450  A4=-l. 45315 

430  A5~=  f.  06  1 41 - 

470  DIM  A(2000) »B(2000) »F(2000) 
430  DTW  ~S'(2000)  F3  (2000) - 


TTT 
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490  FOR  1=1  TO  2000 


510  B  ( I )  =0 


530  FOR  B=1  TO  S 


550  U2=RND(-1) 


B  1  ■  Tl  l  3T4C  lllirtUWl 


570  X3=Z1*C0S<6.233185*U1) +N1 


590  S8=S8+X3 


.•TiflKaroi c«  x : 


610  U3=RND ( -1 ) 


630  Z2=SQR ( -2*D2**2*L0G ( U4 ) ) 


650  IF  Y3<  =0  GO  TO  610 


670  T9=T9+Y3**2 


350  PRINT 


370  D5=SQR( (T3/S)-(X4+*2) ) 


•j 'J v  /iw-g  7  /  g 

390  D6=SQR( (T9/S)- (X5**2) ) 


910  B5= ( 2*N7) **B3 


930  PRINT  "DUCTILITY :AVG»STD  DEV: 


950  PRINT  "APPLIED  STRAIN:AVC»STD  DEV 


,va.  *?-***,. 


i 


* 


PL-1  CCont'd) 


990  PRINT 


■fiT" FIRST  FA ILURE's ,r TINT  fL+ .  5] 


1010  PRINT  “N (MEDIAN) =" 5 INT (N7+ . 5) 


1030  PRINT  "SAMPLE  SIZE="5S 


rniwi 

1050  PRINT  "BETA*" » B1 


150®  riUNI 

1070  D8=SGR(2*B1-(3.14159/Bl)> 
1080  D9= ( (2*N7) **B2) /D8 
1090  P0* <D5**2) + (B5*D6**2) 


Rn 


1110  N9= ( 1/2) * ( (N1/N2)**B1) 


mr 


1130  N=  ( I  - 1 ) +W 

1140  P1=SUR(P0+SGN<  (N/N7)  -  1 )  *09*D5*D6) 
1150  M3= ( X4/P1 ) * ( ( < (N+W) /N7) +B2) - 1 ) 


1170  IF  M3>-0  GO  TO  1190 


iio»  na-*nj 

1190  T2  =  l/ (1+0.2316413#M5) 


1210  C6=(A1*T2)+C5+(A5*T2**5) 


m  W-TII  31^(31  y  LtiyjfMlil 


1230  El  * ( 1 /2) *ABS ( A8) 


IttW  C.^-C.1 

1250  IF  M3>=0  GO  TO  1270 


1270  M4  = (N1/P2) * ( ( ( (N+W) /N9)+B2) -1) 


jL£.«icr  no-nt 

1290  IF  M4>=0  GO  TO  1310 


x  s?vv  no--nt 

1310  T3  =  1/ ( 1+0 . 2316418*M6) 


**3)+(A4*T3**4> 
1330  C3=(A1*T3)+C7+(A5+T3**5) 

13-40  A9*  l  -  (C3*TEXP7 -  ( M5*#2)Y2)TT - 

1350  E3=(1/2)*ABS(A9) 


1370  IF  M4>=0  GO  TO  1390 


1390  F  ( I )  = . 5+E2 


1410  NEXT  I 


1430  FOR  1=2  TO  K 


1450  NEXT  I 

1450  PRINT-"  ">  "TALLIED"  firCffL3rD,,V'"TAn:rED,rT,rCAL(r,D" - 

1470  PRINT  "CYCLES" i "FAILURES" » "FAILURES" » "F (N) " » "F (N) " 


1490  N=  ( I  - 1 ) *W 


xjvyj  1  r  ou;\,773  IfltM  10419 

1510  IF  A ( I ) < .5  THEN  1530 

— rtt^-PRtNT-m-A  1 1  rrlNTCO ( TT+ .  5)  t  S  ( rnTNT(TE4*F ( D  +T5T7  IE 4" 
1530  NEXT  I 

-t540~END - -  — - - 


'’rot 


MP-  t  , 


PLASTIC  REGION  HISTOGRAM  RESULTS 

Table  II  shows  the  desired  and  tallied  parameters  for  eleven  cases  in  the 
low  cycle  fatigue  region.  S  “  2  for  most  structural  materials.  The  sample 
size  of  Nj  for  each  case  is  10,000  to  minimize  the  variances  of  the  results. 
Nm  was  chosen  to  cover  the  upper  and  lower  ends  of  the  low  cycle  fatigue 
regions.  Cases  |  1  j  ,  |  3  |  ,  f  5~|  and  [  7  [  have  d£  *  0.  The  theoretical 

results  for  these  cases  were  rigorously  derived.  The  variances  Ae  and  6g 
were  chosen  in  some  cases  to  be  large  enough  to  cause  the  cycles  to  first 
failure  to  be  significantly  lower  than  Nm.  See  figures  14  -  23. 

The  curves  of  figures  22  and  23  have  the  same  parameters.  The  theoretical 
expressions  for  are  different  however.  There  is  a  better  fit  (especially 
in  the  region  of  first  failures  between  the  tallied  and  theoretical  results 
of  figure  23  where  the  proposed  form  equation  (15)  for  is  used. 

It  can  be  seen  that  there  is  good  agreement  between  the  theoretical  results 
using  the  proposed  fatigue  life  equations  and  the  tallied  Monte  Carlo 
simulation  results. 


31 


TABLE  II 


DESIRED  VERSUS  TALLIED  PARAMETERS: 
LOW  CYCLE  FATIGUE 


■ 

T 

V 

e 

^  c 

m 

DESIRED 

DESIRED 

DESIRED 

DESIRED 

DESIRED 

CASE 

TALLIED 

TALLIED 

TALLIED 

TALLIED 

TALLIED 

S3 

0.2 

0.01 

0.0024 

0 

3472 

0.19976 

0.00996 

1.29  E-5 

3464 

S3 

0.2 

0.01 

0.0024 

0.0003 

3472 

0.19982 

0.01014 

0.002398 

0.0002978 

3471 

S3 

_ 

0.2 

0.03 

0.0024 

0 

3472 

0.19986 

0.0298 

0.0023444 

1.29  E-5 

3467 

0 

0.2 

0.03 

0.0024 

0.0003 

mmsnmm 

0.19948 

0.03029 

0.002348 

0.000298 

^KTTiHHi 

S3 

0.2 

0.01 

0.00632 

0 

0.20008 

0.00999 

0.00632 

4.48  E-5 

i3 

0.2 

0.01 

■EHjlfl&jH 

MM 

0.20005 

0.010003 

»lvlWi 

■B 

m 

0.2 

0.03 

0 

500 

0.199751 

0.03004 

HTimliwttjll 

4.48  E-5 

499 

IE 

0.2 

0.03 

0.00632 

6.00079 

500 

0.20003 

0.03009 

0.006324 

0.000795 

500 

23 

0.2 

0.02 

0.0024 

0.00048 

0.19966 

0.0201 

0.002397 

0.000485 

mSmSmm 

EE 

0.2 

0.02 

0.0024 

0.00048 

0.19964 

0.0203 

0.002395 

0.000478 

HH 

mi 

0.2 

0.04 

0.0024 

0.00048 

3472 

“87X4558 

0.0393 

F  Q.Q02399 

0.000478 

3464 

e,  S^'v-  in/ in 
Nm  'v  CYCLES 
SAMPLE  SIZE:  10,000 
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ELASTIC  REGION  HISTOGRAM  RESULTS 


Table  III  shows  Che  desired  and  tallied  parameters  for  eight  cases  in  the 
high  cycle  fatigue  region.  Three  values  of  3  are  used.  6  *  9.6  represents 
7075-T6  Aluminum  Alloy.  3  »  12.1  represents  G-10  Epoxy  Fiberglass.  S  *  22.37 
represent  AZ31B  Magnesium  Alloy.  Such  a  range  covers  ductile  to  brittle 
materials.  Nm  is  chosen  to  cover  from  the  lower  to  upper  end  of  the  high 
cycle  region.  The  standard  deviations  are  large  enough  to  cause  the  cycles 
to  first  failure  to  be  much  less  than  Nq. 

Figures  25  through  32  show  the  theoretical  results  using  the  proposed  fatigue 
life  equations  and  the  tallied  simulation  results.  Figure  31  is  the  one 
exception. 

I  2  2~^ 

*yA£  +  6£  i3  used  for  case  (T)  ,  not  the  proposed  equation  (15). 
Figure  31  shows  that  the  theoretical  histogram  is  reasonably  accurate  but 
not  nearly  as  accurate  as  those  using  equation  (15)  for  f£.  In  the  vicinity 
of  the  early  failures  the  theoretical  curve  Is  non-conservative.  This  is 
the  reason  that  V*  was  discarded. 

In  all  other  cases  the  proposed  results  are  accurate  compared  to  the  tallied 
results. 

Figures  29  and  32  show  the  effect  of  large  3  on  the  spread  of  cycles  to 
failure.  Figure  32  in  particular  shows  a  preponderance  of  failures  in  the 
early  life  region,  much  less  than 


44 


h 


» 


TABLE  m 


DESIRED  VERSUS  TALLIED  PARAMETERS : 
HIGH  CYCLE  FATIGUE 


£u’  e’  'V  ^  in/in 


N  ^  CYCLES 
m 


SAMPLE  SIZE:  10,000 


45 


CASE 

s 

4  & 

£ 

N 

in 

s 

DESIRED 

DESIRED 

DESIRED 

DESIRED 

DESIRED 

TALLIED 

TALLIED 

TALLIED 

TALLIED 

TALLIED 

© 

0.0185 

0.000925 

0.004387 

0 

500,000 

9.6 

0.0185 

0.000919 
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CYCLES  TO  FIRST  FAILURE  RESULTS 

Appendix  B  shows  the  derivation  of  the  expressions  for 

.  PL-2  is  a  listing  of  the  Monte  Carlo  simulation  program  to  tally  the 

cycles  to  first  failure.  This  program  generates  Nf  samples  in  the  same  way 

as  PL-1.  The  program  first  generates  10,000  samples  and  selects  the  lowest 

valued  sample  as  the  first  sample  of  the  cycles  to  first  failure  N^.  This 

process  is  repeated  19  more  times.  A  total  of  20  x  10,000  ■  200,000  N^  samples 

are  generated  to  obtain  20  samples  of  N^.  and  are  measured  for  these 

20  samples.  The  smallest  single  value  of  N,  is  called  Ni  .  Equations  (17), 

A  ■‘■MIN 

(18)  and  (19)  are  used  to  calculate  the  expected  corresponding  values.  Table  IV 
compares  tallied  and  calculated  results.  Most  of  the  results  show  excellent 
agreement  between  tallied  and  calculated  values.  Some  comparisons  are  good. 

Two  are  poor.  Such  a  range  in  quality  of  agreement  is  considered  to  be  caused 
by  the  low  sample  size  for  N^  of  20  and  not  by  an  inherent  inaccuracy  of  the 
proposed  equations.  To  significantly  increase  the  sample  size  of  would  be 
prohibitive  in  terms  of  computer  time  and  cost.  The  overall  good  agreement 
already  shown  does  not  warrant  any  further  effort. 

PL-2  was  modified  by  changing  lines  130  and  140  from  10,000  and  20  to  400  and 
500  respectively.  Thus  400  x  500  ■  200,000  Ng  samples  were  generated  to  acquire 
500  samples  of  N^.  This  was  done  to  obtain  the  shape  of  the  N^  histogram. 
Figures  33  and  34  show  two  such  histogram  shapes.  The  shapes  look  surprisingly 
like  those  of  only  backwards  (i.e.  rotated  about  the  median  value  of  N) . 


Equations  (17) ,  (18)  and  (19)  are  considered  to  be  accurate. 
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TABLE  I V  COMPARISON  OF  CYCLES  TO 
FIRST  FAILURE  RESULTS 


N 

W1  MIN 

TALL'D 

CALC’D 

TALL'D 

CALC'D 

TALL'D 

CALC'D 

2255 

2300 

2026 

2080 

73 

73 

1398 

1090 

1211 

849 

88 

80 

644 

678 

275 

360 

145 

106 

532 

568 

316 

147 

i 

99 

140 

328 

331 

304 

300 

1 

11 

10 

199 

157 

164 

122  | 

1 

13 

12 

92 

98 

48 

52 

19 

15 

71 

82 

54 

21 

8 

20 

111 

576 

649 

288 

71 

96 

111 

576 

649 

288 

1  71 

96 

158 

69 

13 

9 

i  60 

! 

20 

61,916 

i 

69,356  j 

J 

39,116 

42,749 

8,235 

8,869 

35,177 

32,507  | 

12,374 

i 

6,014 

i  8,807 

8,831 

662,338 

693,556  | 

i 

547,341 

427,421 

84,734 

88,712 

356,562 

326,51l!j  201,893 

60,428 

64,844 

88,694 

615,847 

565,935  i 

l 

410,963 

69,451 

139,375 

165,495 

947,484 

693,326  ! 

800,307 

706,843 

83,542 

-  4,505 

981 

734  ; 

1 

442 

17  i 

1 

389 

239 

All  values  are  in  units  of  CYCLES 
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COMPARISON  WITH  EMPIRICAL  DATA 

The  proposed  fatigue  life  expressions  have  previously  been  shown  to  agree 
well  with  the  Monte  Carlo  simulation  tallied  results.  Now  the  theoretical 
and  tallied  results  will  be  compared  with  the  empirical  results  reported  in 
references  {2},  {8}  and  {10}.  In  reference  {8}J.T.  Broch  describes  fatigue 
life  test  results  of  G-10  fiberglass  single-degree-of-freedom  end  mass 
cantilever  beams  subjected  to  random  stresses.  A  sample  size  of  100  beams 
was  used  for  the  tests.  The  test  parameters  are  as  follows: 


a  -  12.2  ksi;  6  -0.348  ksi;  4  -1.75  ksi 
C  -  33  ksi;  8  -  12.1;  E  -  2700  ksi 
The  corresponding  strain  parameters  are 


-  1/8 

eu  *  2  £  -  0.0129427  In/ In 

E 

£  —  —  -  0.0045  in/ in 

E 


4/  -  0.000647  in/ in 
'E 

<5  ,  «  0.000129  in/ in. 

'E 


Figure  35  shows  a  comparison  of  the  theoretical  and  tallied  histograms 
for  the  above  parameters.  Large  variances  in  the  tallied  are  noted.  However, 
the  overall  shapes  are  in  general  agreement.  Figure  36  compares  theoretical 
and  empirical  data.  Again  large  variances  are  noted  in  the  empirical  data. 

The  overall  shapes  are  in  general  agreement.  Figure  37  compares  the  empirical 
and  tallied  histograms.  They  too  generally  are  in  agreement  with  each 
other.  Figure  38  shows  that  the  variance  of  the  tallied  data  is  smoothed  out 
considerably  as  expected  by  increasing  the  sample  size  from  100  to  10,000. 

This  indicates  that  the  previous  relatively  large  variances  for  the  tallied 
and  empirical  data  is  an  expected  result  of  the  small  sample  size  of  100. 
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Figure  39  is  a  histogram  of  N^.  Also  shown  are  calculated  and  tallied 
values  of  and  empirical.  Quantitatively 

*  33,983  cycles 

cal'd 

Ni  ■  37,914  cycles 

Tallied 

»  42,950  cycles 

empirical 

All  of  the  data  indicates  good  agreement  between  theoretical,  Monte  Carlo 
and  empirical  results. 

Figure  40  shows  additional  empirical  fatigue  failure  data  {10}.  Again  the 
theoretical  results  are  in  good  agreement  with  empirical  results. 
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SYMBOLS 


A  material  constant;  true  ultimate  stress 

b  fatigue  strength  exponent 

c  fatigue  ductility  exponent 

C  constant  of  random  fatigue  curve 

E  modulus  of  elasticity 

erf (a)  error  function  of  argument  a 

F(N)  probability  of  failure  at  N  cycles 

N  applied  stress  cycles 

Na  histogram  bin  width 

9 

Nf  number  of  stress  cycles  to  failure 

%  median  stress  cycles  to  failure;  cycles  to  50Z  failures 

stress  cycles  to  first  failure 
average  value  of 
minimum  value  of 
N8  random  variable 

p(a)  probability  density  function  of  a 

q  histogram  quantity 

rms  root  mean  square 

S  total  sample  size 

AS  applied  sinusoidal  "engineering"  stress  amplitude 

2 

x,y,z  general  variables 

r,h,v  i 

a 


3 

r(a) 


general  variable 

fatigue  curve  slope  parameter 

gamma  function  of  argument  a 


A 


standard  deviation  of  stress  fatigue  curve 
standard  deviation  of  strain  fatigue  curve 
standard  deviation  of  applied  stress 
standard  deviation  of  applied  strain 

applied  strain  amplitude  (one-half  applied  strain  range) 
ultimate  strain  amplitude;  ductility 


average  value  of  e 


£U 

4 

&£ 

A£< 

Ae 


Aa 

~T 


average  value  of 

fatigue  ductility  coefficient 

applied  strain  range 

applied  elastic  strain  range 

applied  plastic  strain  range 

correction  factor 

fatigue  strength  coefficient 

average  value  of  random  rms  stress 

applied  sinusoidal  "true"  stress  amplitude 

standard  deviation  of  N, 


resultant  stress  standard  deviation 
resultant  strain  standard  deivation 
modified  form  of 
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APPENDIX  A 


DERIVATION  OF 


p(Nf> 


From  equation  (1) 


(A-l) 


This  is  of  the  form  z  *  y 


where  x  *  e, 


y  ■  ek 


z  -  2N£ 


x  and  y  are  simultaneous  random  variables 


Let  £p  be  a  Gaussian  random  variable  with  average  ey  and  standard  deviation 


A£.  From  reference  {9} 


f(eu)  -  x  -  eu 


eu  -  x 


d 

d  e„ 


P(eu) 


(A-2) 
(A— 3) 


(A-4) 


(A- 5) 


JM 

iJ 


% 

i 

% 


1/8-1 


8 


Similarly 


1/6-1 


P(y)  -  - 


)  l/B  _ 
!X  -  eu 


8  5e  2tt 


.  (  v1/8-  gV 


2  Je  - 


(A-6) 


CA-7) 
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APPENDIX  A  (Coat’d) 


P„  „  Cx.y)  -  p(x)  p(y) 
$»7 


Jfe.y  Cx»y) 


p(z)  -  2 


1/B-l  1/ 0—1 


8  A£  5e  2* 


y  p_  „  (z  y.  y)  dy 


(  1/8  _  2  1/8  _  2 
1  J  (x  -  £u)  +  (y  -  £u) 


After  much  manipulation  It  can  be  shown  that 


1/8  1/8-1  -  ft . IT- v> 

P(Nf}  ‘  8  AeN!‘“  ”  2  VF  e  r  f  2  “  «£(«!> 


i  “hVr  /2  hA  . 

_e  .( -  erf  (a2) 


1  (2Nf) 

where  r  ■  "J  -  ^ 


+  t2 


»•-  -ia>-  *  f2 

Ae  e 


1  e 


V  “  2 


Ctx  -  a/T  f  £  +^p 

IT 

a_  ■  h//r 


(A-8) 


(A-9) 


(A-10) 


APPENDIX  A  (Cont'd) 


0 


erf  (0)  ■  0  ;  erf  (“)  *0.5 

erf  (-a)  ■  -  erf  (a) 

Equation  (A-10)  Is  the  expression  for  p(N^)  when  and  e  are  simultaneous 
random  variables. 
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APPENDIX  B 


DERIVATION  OF  Nx 

N^  -  average  value  of  cycles  Co  firsc  failure 
S  -  number  of  opportunities  for  failure 
S  -  10,000  for  these  cases 


F(NX)  -  0.5  +  erf  (a);  erf  (a)  -  -0.4999 


a  -  erf  (-0.4999) 


-3.7195451 

-3.7195451 


Solving  for  Nj. 


N1  *  Nm 


-T  s 


1  - 


3.7195451 


N. 


1MIN 


1  i  eu  -  4.52As 
7  +  4.526. 


%  — 


%  -  (Nx  -  Nx  )/3 


MIN 


(B-l) 


(B-2) 


(B-3) 


(B-4) 


(B-5) 
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